STABILITY OF SPECTRAL EIGENSPACES IN NONLINEAR 
SCHRODINGER EQUATIONS 
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' Abstract. We consider the time-dependent non linear Schrodinger equations 

with a double well potential in dimensions d = 1 and d = 2. We prove, in 
the semiclassical limit, that the finite dimensional eigenspace associated to the 

' lowest two eigenvalues of the linear operator is almost invariant for any time. 

< 

CO . 

1. Introduction 

Here we consider the time-dependent nonlinear Schrodinger equation (hereafter 
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NLS) 



dt ' 
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. where 

it (2) h °=-£i a+v > A =E^ d=i ' 2 ' 2 ™ =i 

i-^h ' is the linear Hamiltonian operator with symmetric double- well potential V(x), and 

d • where 

S : (3) W = W{\i>\ 2 ) = \4,\ 2a , *>0 

is a nonlinear perturbation. 

We recall that a double well potential is a positive potential symmetric with re- 
jJJ ■ spect to the reflection through a hyperplane, and having two nondegenerate distinct 

absolute minima. When the nonlinear term is absent, the linear Hamiltonian Hq 
has the two eigenfunctions corresponding to the two lowest eigenvalues which have 
either even or odd-parity. Given an initial datum which is a linear combination 
of such lowest eigenfunctions, then the corresponding solution performs a beating 
motion, namely the probability density oscillates periodically from a state almost 
concentrated on one minimum to a state almost concentrated on the other one. 
The beating period usually plays the role of unit of time. When the nonlinear term 
is restored, a symmetry breaking phenomenon occurs: that is, if the strength e of 
the nonlinear term is larger than a threshold value then new asymmetric stationary 
states appears pQ, 0], Furthermore, for higher strength e of the nonlinear 

term, the beating motion is generically forbidden [5], 0, HD- These results can be 
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hcuristically obtained by reducing the NLS equation to a 2-dimensional dynami- 
cal system, namely the restriction of (QJ to the space generated by the lowest two 
eigenvectors. Mathematically the difficulty consists in proving that the dynamics 
of the complete NLS is close to the dynamics of the 2-dimensional reduced system. 
In §|, making use of semiclassical estimates and refined existence results for NLS 
this stability result has been obtained for times of the order of the beating period 
in dimensions d — 1, 2 and any a £ M + . 

In the present paper we concentrate on the case of local nonlinearity © with a 
a positive integer, in this framework we extend the previous result by jHj proving 
that, in the semiclassical limit, the 2-dimensional eigenspace is almost invariant 
for any time (Theorem However, due to the possible presence of positive Lya- 
punof exponents, our result does not allow to show that the 2-dimensional system 
describes the dynamics over time scales larger than e _1 . 

Our result is obtained here making use of variational methods and by introducing 
the scale of Hilbert spaces X s = T>{H^ 2 ), s > 0, constituted by the domains of the 
powers of Hq, endowed by the graph norm; in fact, in dimension d — 1 and d = 2 
we make use of the energy space X\. 

We close this section by introducing some notations: 

- The notation y = 0(e~ r / h ) means that there exist h* > and a positive 
constant C > 0, independent of H, such that 

\y\ < Ce'V*, Vfi £ (0, H* ) 

- The notation y — 0(e~ r / h ) means that for any L', < V < L, then y = 
0(e~ r / h )] that is, there exist h* > and a positive constant C — Cj~> > 0, 
independent of h, such that 

\y\<Ce- T '/\ Vh€(0,H*). 

- As usual, || • || lp usually denotes the norm of the space L p (M. d ): \\-\\ s denotes 
the norm of the Hilbert spaces X s , in such a notation Xq = L 2 (M. d ) and 

II • llo = II • IU- 

- As usual, R denotes the set of real numbers, M + denotes the set of positive 
real numbers, N denotes the set of positive integer numbers, C denotes 
any positive constant independent of h and t, and (x) = yl + |x| 2 , |a;| 2 = 
52j=i x "j: where x = (x%, . . . ,x d ) £ R d . 

2. The model: preliminary assumptions and main result 
Hereafter we assume the dimension d — 1 or d = 2. 

2.1. Double well potentials. 

Hypothesis 1. The double-well potential V(x) £ C°°(R d ) is a symmetric real 
valued function such that: 

i. V(x) admits two minima at x — x±, X- ^ x + , such that 
V{x) > V mm = V(x±), V.t £ R d , x + x ± ; 
for the sake of definitess let us assume 



mitt 



= 1. 
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ii. There exists m > such that, for large \x\ one has 
C(x) m < V(x) < C-\x) m 



for some constant < C < 1; 
iii. For any multi-index a £ N d there exists a positive constant C a such that 

\d^d^---d^ d V(x)\<C a (x) m -^ 

Remark 1. For the sake of definiteness, we assume that the symmetric potential 
is such that 

V(-X1,X2, ■ ■ ■ ,Xd) — V(X1,X2, ...,Xd) ■ 

Furthermore, for the sake of simplicity we assume also that 

WV(x ± ) = and Hess V(x±) > 0; 

the case of degenerate minima, that is det[HessV (x±)\ = 0, could be treated in a 
similar way; however, we don't develop here the details. 

Remark 2. As an example one can consider the following double-well potential 



2.2. Well-posedness of the Cauchy problem. Under Hyp. 1 the Cauchy prob- 
lem is locally well posed. 

In fact, in the quadratic (i.e. m = 2) and sub-quadratic (i.e. m < 2) cases if 
V>° S H 1 is such that 



then (see JT] and §9.2 in [3]) the Cauchy problem admits a unique bounded local 
solution -0* £ H 1 such that 



d 



V(x) = 1 + ^2^x1 + Vi(x) 



where V\ £ -L°°(R d ) is a symmetric function 



V 1 (-x 1 ,x 2 , ...,x d ) = V 1 (x 1 ,x 2 , ■ ■ - ,x d ). 





for any t £ [0, 6), for some S > 0. 

Moreover, ijj 1 satisfies the following conservation norm 

(4) M{^)=M{^), 



where 




and 
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where 

foWH^W)^ and P(^) = -^||V^ +1 ||2 2 . 

1 + cr 

For super-quadratic potentials (i.e. m > 2) if ip° G X s then (see Theorem 1.5 in 
the Cauchy problem admits an unique solution ip* for any t G [0, S), where 5 
continously depends on ||?/> || s . Furthermore, the solution ip* continuously depends 
on t and on the initial condition ip° . The conservation of the norm and of the 
energy JSJ) will follow from standard density arguments 3 and from the Gagliardo- 
Nirenberg inequality, which implies, in dimension d = 1 and d = 2, the continuity, 
with respect the norm H , of the energy functional: 

Notice that when e > then the conservation of the energy implies that || V^|| < 
C for any time, hence the solution is global, i.e. 6 = +oo. For e < 0, in the 
sub-critical case no < 2 then the conservation of the energy still give the same 
estimate on ||V-0|| and the global existence of the solution. In fact, for e < small 
enough, by energy conservation and bootstrap arguments (see any solution is 
global in the critical no = 2 and hyper-critical no > 2 cases too. 

2.3. Linear beats. The operator Ho formally defined by J2J admits a self-adjoint 
realization (still denoted by H ) on L 2 (R d ) (Theorem EEL 1.1 in with purely 
discrete spectrum o(Hq) = o c i{Hq), where Od(Ho) denotes the discrete spectrum. 
Let Afc, k E N, be the eigenvalues of Hq 

Ai < A 2 < A 3 < A 4 < . . . < A fe < . . . 

with associated normalized (in L 2 ) eigenvectors <fik[ x )- I n particular, the lowest 
two eigenvalues of Hq are non-degenerate and there exists < C < 1, independent 
of H, such that 

1 + Ch < Ai 2 < 1 + C- X h and inf [A - Ai 2 ] > CK 

Ae<r(H a )-{A 1 , A 2 } 

and 

(6) A - A 2 > CXti, VA g o(H ) - {\ u A 2 } . 

It is well known that the splitting between the two lowest eigenvalues 

(7) w = i(A 2 -A0 

vanishes as H goes to zero. In order to give a precise estimate of the splitting 10 let 
T = inf f ^V{x)- V mm dx > 

be the Agmon distance between the two wells; where 7 is any path connecting the 
two wells, that is 7 G AC([0, l],K d ) such that 7(0) = X- and 7(1) = x + . From 
standard WKB arguments (see 0) it follows that the splitting is exponentially 
small, precisely that 

(8) uj = 0(e- r / h ). 

The normalized eigenvectors (^1,2 associated to Ai, 2 can be chosen to be real- 
valued functions such that ipi and <p 2 are respectively of even and odd-parity: 

<Pj{-xi,X2, ...,x d ) = (-iy +1 (pj(xi,X2, ■ ■ -,x d ), j = 1,2; 



5 



Wc define now the single well states 

PR = [<Pi + <P2] and <p L = -j= [ifx - ip 2 ] 

such that 

<Pr{-X\,X2, ...,Xd) = Vl(xi,X2, ■■■,Xd) 

They are localized on one well in the sense that for any r > there exists C r > 
such that 



/ \v R (x)\ 2 dx = l+0(e- c "/ h ) 

JD r (x+) 



and 



\<PL{x)\ 2 dx = 1 + 0(e-°r/ h ) 

£>,(*-) 

where D r (x±) is the ball with center x± and radius r. For such a reason we call 
them single-well (normalized) states. In particular 

\\<Pr<Pl\\l- = 0{e- T ' h ) 

Let 

n = (pi, -)l 21 Pi + (<P2, -)l 21 P2 and IL C = I - II 

be the projection operator onto the eigenspace orthogonal to the bi-dimcnsional 
space associated to the doublet {Ai^}- We will study the dynamics of equation Q 
with initial data almost exactly concentrated on <I>o :=Span(^i, if 2)- 

Remark 3. If one takes initial data ?/> in $0, i-e. of the form 

4> = C1V1 + C2P2 = CrPr + ClPl, Cr = ^= , Cl = ^= 
then the linear Schrodinger equation 
ihi) 1 =H tp t , 

has an explicit solution given by 

ip*(x) = e ' int/h [(CrPR + Clpl) cos(tut/h) + i (Clpr + (r<pl) sm(wt/H)] 
where 

/n\ Ai + A 2 A 2 - Ai 

(9) n = ^—> "=^— 

That is ^)'(x) performs a beating motion with beating period 

T = 

Such a period usually plays the role of unit of time. 
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\nL^<—\\(- h2 ^) 1/2 n%-u\\^ e ° 



2.4. The nonlinear system. To introduce the analytic framework in which we 
will work we first give the following 

Definition 1. For any integer number s > we define the Hilbert space X s := 
D(Hq 2 ) equipped with the norm 

\m 2 s -.= \\Hfni^m^) L * 

Remark 4. Here we will use the energy space X\ . In particular, the Gagliardo- 
Niremberg inequality yields the following 

C_ 

where \\^p\\ L 2 = 1 and e a a suitable positive number. Hence, we can conclude that 

W\\». + * < 

since 

(-H 2 A) h iP 2 = (yj,-h 2 A^) < (^,H^) = || 

L 2 

3. Main Results 

Hypothesis 2. Let u) be the splitting (JJJ (which satisfies the asymptotic estimate 
|5| ), and let e be the strength of the non-linear term. We assume that the real- 
valued parameter e depends on h in such a way 

\ Aft- da/ 2 

(to) ^ <c, Vhe(o,h*) 

to 

for some positive constant C, independent of H, and for some H* . 
Remark 5. The ratio 

e ft-da/2 



V 



to 

plays the role of effective nonlinearity parameter. The above assumption implies 
that \r)\ < C. 

The main result of this section is the following 

Theorem 1. Assume Hypotheses 1,2, and consider the Cauchy problem (|T|l -(|3*)l . 

Then there exist positive constants < C < 1 and 7 such that, if H is small enough 
and the initial datum ip fulfills 

(11) l|n c ^>°|| 1 < ciel 1 / 2 , 

then one has 

(12) l|n c v/||i<^, vteR. 

Hence the 2-dimensional space $0 := n [L 2 (M d )] is almost invariant. Thus 
one expects that corresponding to initial data satisfying l|llf) the dynamics is well 
described by the restriction of the equations of motion to $0. Actually such a re- 
stricted dynamical system coincides, up to (formal) order e 2 , with the 2-dimensional 
dynamical systems 

, > ( ihc R = -ujc l + flc R ; + eC a \c R \ 2a c R 

[ ' 1 ihc L = nc L -u;c R + eC a \c L \ 2 °c L ' 
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where 

c* = \\<(%\\h = = o(h- d °' 2 ), 

which has the integral of motion 

l(c R , c L ) = n(\c R \ 2 + \c L \ 2 ) - lo (c R c L + c L c R ) + C a -^— (|c fl | 2(CT+1) + |c L | 2(ff+1) ) 

a + 1 V / 

Moreover, such an Hamiltonian system has an independent integral of motion (the 
quadratic part of the Hamiltonian), and thus it is integrable. The 2-dimensional 
system was studied in detail in 5 obtaining that when the nonlinearity parame- 
ter r\ = ^—^ — - is large enough almost all its solution do not posses the beating 
property, i.e. the probability density remains concentrated in one well. 

Concerning the relation between the solution of the two dimensional system and 
the solution of the complete system we have the following 

Corollary 1. Let f a (t) — C R ip R + Cl^l, where C R and Cl are the solution 
of the 2-dimensional system with initial datum <p a (0) with of norm 1, i.e. 
\ c r\ 2 + \ c l\ 2 = 1- Denote tp(t) := Hip(t), where ip(t) is the solution of the complete 
system (JTy) with the same initial datum. Then there exist a positive C and a positive 
a such that the following estimate holds: 

y a (t)-^t)\\ L2 <c e ^\t\ . 

4. Proof of the main results 

Let 

OO 

k=l 

and define the Hilbert spaces £ 2 of the complex sequences £ = {Cm} m gN sucn that 

iicn^E A ^-i 2 <^ 

and remark that in such a way we have defined the correspondence 

V> G X. «-» C = U(1>) G l\ 

which is a unitary isomorphism. 

In terms of these variables the norm and quadratic part £ o of the Hamiltonian 
are given by 

Af = ]TKfc| 2 = IICIIo and £ = Y,^\Ck\ 2 HICIIi- 

k>l k>l 

Remark 6. In order to simplify the notations, we rescale time by the trans- 
formation tt-^t/h. With such a notation the beating period is now given by 

T = — 

UJ 
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4.1. Variational results. We recall that we work here in the energy space X\. 
Thus, in this section, norms and distances will always be in this space. 

To exploit the fact that both £ e and M are conserved along the flow we have to 
introduce a few geometrical lemma. 

Denote 



wow := ttoai 2 + IC2I 2 ) + £ A fc iai 2 = £ 1*10 

fe>3 fe>l 

and 

:=^(|C 2 | 2 -|Ci| 2 ) + ^W , 

where v\ — V2 = fl, Vk = Afc for k > 3 and fl and w are defined in JJ]J, so that 
£ e = H + V e . 

In this section we will almost always use real coordinates pj , qj defined by 

_ qj + ipj 
* s/2 ' 

so all the functions will be considered as functions of tp, (p,q), or Q according to 
convenience. 

Define the (smooth) surface 

S := {ip G X x : N{4>) = 1} , 

and consider the function ho := Wo| 5 . 

Lemma 1. The manifold 

N :=SnU(L 2 (R d )) = 

ip G Xi : ijj = Cifi + (2^2 , \Ci\ 2 + \C2\ 2 = Pl+ql+ / 2 + ql =1 

is an absolute minimum of ho. Furthermore, for any point ip € N decompose the 
tangent space T^S as 

T^S = T. 4 ,N®(Tt P N) ± . 
Then, the second differential (&ho of ho at a point ip G N is such that 

(14) d%ho{x,x) > Ch\\x\\j , vx e (t^n) 1 - 

Remark 7. Since the function ho is smooth one also has 

(15) &%h {X,X) < C\\X\\l , VX G (T^N) 1 - 
with a suitable C . 

Proof. It is a trivial application of the method of the Lagrange multipliers. Con- 
sider Hq + XAf; the critical points of ho are obtained by finding the zeros (on <S) of 
the equations 

fin d(H + AAQ I 

(16) ^ = 2^+A> j) 

where Vj = Q for j = 1, 2 and i/j = Xj for j > 3, and of the analogous equation for 
<7j. Thus a solution of l|16[l is given by A = — f2 and tp E N. It follows that TV is 
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constituted by critical points of ho. It follows that the second differential of ho at 
such points is well defined. Moreover, we recall that, given a vector 

Y = ^2(Q j + iP j )i Pj eT^S, 

one has 

+ 2P)0 ,2M + w ,^p) 



3,1 



= $> + A)(P 2 



j>i 



where A = — Q is the Lagrange multiplier determined from the criticality condition. 
By the condition Y € (T^N) 1 - one has Pj = Q 3 ■ = for j = 1,2 and thus one has 

d 2 /io(^, F) > ^ Cft (P? + Q 2 ) = Ch\\Y\\\ , vr e {T^pN) 1 - 
i>3 

since JSJl, and therefore the thesis on the second differential follows. In particular 
N is a minimum of h^. In a similar way one can show that all the other critical 
points of ho are saddle points. □ 

Remark 8. By definition it follows that 

h \ N = n. 

Lemma 2. There exists a positive C such that, provided 

(17) h o (iP)-n<±h , for some [3 > 3, 
£/ien one /ias 

(18) - [d(V, AOf < h (i>) -n<C [dty, AT)f • 

Proof. This is a standard result in differential geometry; we give here its simple 
proof for the sake of completeness. Actually it is based on the use of the exponential 
coordinates (see appendix \K\ for their construction), which are coordinates 

1&-> (n,w) eNx (T^N) 1 - 

in which the intersection of N with the domain of definition of the coordinates 
coincides with w = 0; moreover for any ip = (n, w) one has 

d(tp,N) :=d((n,w);N) = \\w\\i . 

where the distance is as usual defined as the length of the shortest geodesic from ip 
to N. Using these coordinates, consider the Taylor expansion in w, at a point of 
N, of the function 

(19) h (n,w) - h Q (n,0) = d 2 {n>0) h {w , w) + 0(\\w\\l) , 

where ho(n,w) = ho(tp) and ho(n,0) = Q = ho\ N . From this fact and from 
(|f 4J and t| 1 7JI then the left hand side inequality (|f 8(1 follows. The right hand side 
inequality of (fT£|l follow from l(T5|) and lfT^|l . □ 
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Corollary 2. Let tp £ S such that ||II c ^>||i = fj, = 0(h@/ 2 ) for some [3 > 3, then 
(20) Ch?n<d{ip,N) <Cfrifi 

for some C > 0. 
Proof. Indeed, since 

hoW = n(|Ci| 2 + IC 2 | 2 ) + 5> fe |c fe | 2 



fe>3 



then 



= fMICil 2 + IC 2 | 2 ) + l|n c ^ 
= f!-ft||n c ^|| 2 + ||n c ?A|| 2 

&„(v) - n = l|n c vil 2 - nijiie^llg < ||n c v|| 2 = m 2 



From this fact and from Lemma [5] the right hand side inequality ()20[1 follows. 
Similarly 

hoty) - n = ^(Aj - n)|0l 2 > ct^a^I 2 > c^||n c ^l!f 

3>3 i>3 

from which the left hand side inequality (1201) follows. □ 

Proof of theorem ^ First remark that ip° € iS and 1 1 n c '0 1 1 1 — M — 0(e 1 / 2 ) 
implies (see Corollary^ d{iP°,N) < Chr 1 ! 2 ^ = Oie 1 ' 2 ). Then, remark that from 
GDI then 



£«> ) = ||^ ||2 + |j^+l||2 2 



e 

||n^o|| 2 +0(e 1 / 2 ) = fi+0(e 1 / 2 ) 



and that, by conservation of energy and Remark 0] one has ||V>(0lli < 2S7 for all 
times. From this fact and from Remark 0] one also has the a priori estimate which 
holds for dimension d = 1 and d = 2 

IW(*))I < c-^ 

for some /3; thus, using l|18|) one has the chain of inequalities 

[d 2 (m,N)} 2 <jr[hoMt))-n] 

< j [MO - fi) + (s e (m) - £^ )) + - PcfsKt)))] 

< j [(h (i; ) - n) + - P e (V(*)))] < ■ 

From this fact and from i|20|) then the thesis follows. □ 
Finally, Corollary ^ immediately follows by comparing the two-level approxima- 
tion (|13fl with GJ and by means of estimate H12|) and standard Gronwall's Lemma 
arguments (see e.g. [S])- 
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Appendix A. Construction of exponential coordinates 

Let M and N be Riemannian manifolds modeled on Hilbert spaces H and K, 
with norms || • \\h and || • ||#; denote by g the metric of M; let i : N — > M be a 
smooth isometric embedding. In the following, for the sake of simplicity, we will 
identify N and with i(N), and similarly for related objects and spaces (as TN). 

Lemma 3. Let xq be a point of N and let W = (T^^N)- 1 , then there exists a 
coordinate system of M defined in a neighborhood 14 C M of Xo with the following 
properties: 

i. U 3 x h-> (n,w) € N x W; 

ii. let d(x, N) be the distance usually defined as the length of the shortest geo- 
desic from x € N to N ; then d((n, w); N) = \\w\\k- 

Proof. We proceed in some steps. To start with we choose a coordinate system 
(n, wi) with origin ij e JV such that, if (n,wi) £ Vi C N x W, where Vi is a 
neighborhood of xo ; then the intersection of N with the domain of definition of the 
coordinate system coincides with the set (n, 0). Let 

n„:ViCiVx (T {n>0) N) X 

be the orthogonal projector with respect to the scalar product 5( nj o)- Define the 
map 

Vi 3 (n,wi) h-> (n,H n Wi) G V . 

If Vi is small enough, then such a map is an isomorphism on its image, since its 
differential at the origin is invertible. Thus it defines a new coordinate system. 
Let (n, 0) G Vi, then in these coordinates one has 

T (nt0) N = {(X 1 ,X 2 )e NxW : X 2 = 0} 
(T (n , 0) iV) x = {(Xi,X 2 ) eNxW : X x = 0} 

We use such a coordinate system in order to define the needed coordinate system. 
Take (n,w) € N x W small enough, and consider the geodesic "f( n ,w) (s) starting 
from (n, 0) with initial velocity w. Consider the (exponential) map 

(n,w) i-> 7( n , w )(l) , 

by implicit function theorem it is locally invertible (its differential at the origin is 
the identity), and thus any point a: of a neighborhood hi of £o can be represent 
uniquely by the points (n, w) such that 7(„ jt0 )(l) = £o- Remark that moreover one 
has 

(21) t{l)-= I Vs(7(s),7(s))^ = Hk 

Jo 

since 7(0) = \\w\\k and g{ A /(s), A /(s)) is independent of s along geodesic lines. 

In the coordinate system just introduced, let x = (n, w) then its distance from 
N is the length of the shortest geodesic joining x to N. In turn such a geodesic is 
perpendicular to N. Thus (n, 0) is the point where it starts and w is the tangent 
vector at such a point. Then property ii. follows from l|21|) . □ 
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